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Entropy in nonequilibrium statistical mechanics is investigated theoretically so as to extend the 
well-established equilibrium framework to open nonequilibrium systems. We first derive a micro- 
scopic expression of nonequilibrium entropy for an assembly of identical bosons/fermions interacting 
via a two-body potential. This is performed by starting from the Dyson equation on the Keldysh 
contour and following closely the procedure of Ivanov, Knoll and Voskresensky [Nucl. Phys. A 672 
(2000) 313]. The obtained expression is identical in form with an exact expression of equilibrium 
entropy and obeys an equation of motion which satisfies the //-theorem in a limiting case. Thus, 
entropy can be defined unambiguously in nonequilibrium systems so as to embrace equilibrium sta- 
tistical mechanics. This expression, however, differs from the one obtained by Ivanov et at, and 
we show explicitly that their "memory corrections" are not necessary. Based on our expression of 
nonequilibrium entropy, we then propose the following principle of maximum entropy for nonequilib- 
rium steady states: "The state which is realized most probably among possible steady states without 
time evolution is the one that makes entropy maximum as a function of mechanical variables, such 
as the total particle number, energy, momentum, energy flux, etc." During the course of the study, 
we also develop a compact real-time perturbation expansion in terms of the matrix Keldysh Green's 
function. 



I. INTRODUCTION 

Much effort has been directed towards extending equi- 
librium thermodynamics and statistical mechanics to 
open nonequilibrium systems with flows of particles, mo- 
mentum and/or energyjii^ii^iii Beyond the linear-response 
theory,"SjS however, there seems to have been yet no es- 
tablished theoretical framework comparable to the equi- 
librium one^IiS The purpose of the present paper is to 
make a contribution to this fundamental issue, especially 
on steady states without time evolution, by studying the 
roles of entropy in nonequilibrium statistical mechanics. 

The present approach with entropy is motivated by a 
couple of following observations. First, equilibrium sta- 
tistical mechanics is constructed on the principle of equal 
a priori probabilities, and the equilibrium is identified as 
the one which is most probable. It is hard to imagine that 
this concept of "maximum probability" loses its validity 
as soon as the system is driven from outside. For exam- 
ple, a gas initially prepared in one half of the container is 
expected to expand over the whole available space even 
in the presence of heat conduction through the surface. 
In equilibrium, it is entropy that embodies "maximum 
probability," whose statistical mechanical expression is 
given by Boltzmann's principle)^ 



S'oq = fee log W . 



(1) 



And all the other free energies stem from through the 
mathematical procedure of Legendre transformations, 
thereby inheriting the extremum property of entropy. 
Thus, one may expect to have an appropriate description 
of open steady states by extending the concept of entropy 
or "maximum probability" to nonequilibrium situations. 
Note in this context that eq. H]), which represents the 
principle of equal a priori probabilities, is essentially an 
equilibrium expression with no dynamical equation at- 



tached to it. 

The second motivation originates from the classical 
Boltzmann equationi^ i^°i^^i^^ Let us define the total en- 
tropy of dilute gases by 



(2) 



where / — f{p, r, t) denotes the distribution function.— 
Following the procedure to prove the iJ-theorem^i^ we 
then obtain the inequality: 



d»5' f 

I V- js(r-)dV>0, 



(3) 



-^^/(log/ 



1). It hence fol- 



with js{r) = -fee / (2^-. 

lows that dS/dt > for the isolated system. This is 
the usual 7J-theorem. Looking at eq. ^ more carefully, 
however, one may notice that d5/dt>0 holds whenever 
J 'V ■js{r)d^r = 0. Thus, entropy is expected to increase 
monotonically even in open systems as long as there is 
no net inflow or outflow of entropy through the bound- 
ary, besides those of energy, momentum and particles. 
This observation suggests that the open steady state, if 
it exists, may also correspond to the maximum of entropy 
with appropriately chosen independent variables. In this 
context, eq. ([2]) is superior to eq. ([T]) in that it is appli- 
cable to nonequilibrium systems, but inferior to eq. ([T]) 
in that it is good only for dilute classical gases. 

Now, the main purposes of the present paper are 
twofold. First, we derive an expression of nonequilibrium 
entropy for an assembly of identical bosons/fermions in- 
teracting via a two-body potential so as to be com- 
patible with equilibrium statistical mechanics. Such 
an investigation was performed recently for the contact 
interaction in a seminal paper by Ivanov, Knoll and 
VoskresenskyJ^ii^ We here extend their consideration to 
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general two-body interactions, critically reexamine their 
derivationfi^ and present an expression of nonequilib- 
rium entropy which differs from theirs in an essential 
point. To be more specific, we adopt the nonequilib- 
rium Dyson equation for the Keldysh matrix^^ as our 
starting point, which is transformed into a tractable 
form by the gradient expansion, i.e., the procedure well- 
known for a microscopic derivation of quantum trans- 
port equations fi^ii^ii^iiiii^ii^i2£i2ii2^i2^ An expression of 
noncquilibrium entropy density is then obtained from the 
reduced Dyson equation as eq. below. It is a di- 
rect extension of eq. ^ to include both quantum and 
many-body effects in nonequilibrium situations, which is 
also compatible with the equilibrium expression.— We 
will show explicitly that "memory corrections" of Ivanov 
et al.,^^ which is the origin of the above mentioned dif- 
ference in nonequilibrium entropy, are both unnecessary 
and incompatible with equilibrium statistical mechanics. 

Second, we propose a principle of maximum entropy 
for nonequilibrium steady states in §IV CI A key point 
is that we choose mechanical variables, such as the total 
particle number, energy, momentum, energy flux, etc., as 
independent variables of entropy. Indeed, temperature, 
pressure and chemical potential, which are not adopted 
here, are all equilibrium thermodynamic variables defined 
with partial derivatives of eq. ([!]); thus, they cannot spec- 
ify any nonequilibrium state of the system. The principle 
may enable microscopic treatments of open steady states 
in exactly the same way as equilibrium systems. Its va- 
lidity can only be checked by its consistency with exper- 
iments, as is the case for the principle of equal a priori 
probabilities in equilibrium statistical mechanics. Thus, 
it will be tested in the next paper on Rayleigh-Benard 
convectio n'^i^^i^^'^''i^^'^^i'^° of a dilute classical gas, which 
may be regarded as the canonical system of nonequilib- 
rium steady states with pattern formation. It may be 
worth emphasizing at this stage that the present princi- 
ple is connected with entropy itself. Thus, it has to be 
distinguished from the principle of excess entropy pro- 
duction by Gransdorff and Prigogincp which has been 
criticized by Graham,'^^ for example; see also ref. fs^ 

During the course of study, we also develop a com- 
pact perturbation expansion on the Keldysh contourii^ 
In principle, this expansion can be carried out for the 
round-trip Keldysh contour in the same way as in the 
equilibrium theoryi^i^ii^ When writing it with respect 
to the real-time contour of — oo < i < oo, however, one 
usually has to introduce additional contour indice a^^i^^ 
which make the actual calculations rather cumbersome 
and complicated. The present method will enable us to 
carry out the expansion on the real-time contour directly 
in terms of the 2x2 Keldysh Green's function without 
using the contour indices. Among various approxima- 
tions in the perturbation expansion, we here specifically 
consider Baym's ^-derivative approximation.— It has at 
least the following advantages: (i) it includes the exact 
theory; (ii) various conservation laws are automatically 
obeyed; (iii) the vertex corrections, or the Landau Fermi 



liquid corrections in a different terminology, are natu- 
rally included; (iv) n-particle (n = 2,3, •••) correlations 
can also be calculated within the same approximation 
scheme, i.e., there is a definite prescription here to treat 
the Bogoliubov-Born-Green-Kirkwood-Yvons (BBGKY) 
hierarchyii^ The derivation of nonequilibrium entropy by 
Ivanov et ali^ will be reexamined critically within the 
present expansion scheme. 

This paper is organized as follows. In §2, we develop 
a compact real-time perturbation expansion in terms of 
the matrix Keldysh Green's function for an assembly of 
identical bosons/fermions interacting via a two-body po- 
tential. We consider the ^-derivative approximation in 
detail to write down the Dyson equation for the Green's 
function and the expression for the two-particle corre- 
lation function. In §3, we first introduce the spectral 
function A and the distribution function (/> in the Wigner 
representation; they form alternative two independent 
components of the Keldysh Green's function. We then 
carry out the first-order gradient expansion to the Dyson 
equation to obtain the equations for A and 0. In §4, 
we derive an expression of noncquilibrium entropy as eq. 
(IM)) below. A detailed discussion will be given on the 
difference between the present expression and the one 
obtained by Ivanov et al}^ We then propose in ^IW CI a 
principle of maximum entropy for nonequilibrium steady 
states. Section 5 summarizes the paper. In AppendijXl 
we show with the present perturbation-expansion scheme 
that various conservation laws are automatically obeyed 
in the ^-derivative approximation. AppendijBl presents 
expressions of the vertex functions in the second-order 
derivative approximation. In AppendijCl we derive basic 
conservation laws in the first-order gradient expansion of 
the ^-derivative approximation. Finally in Appendi>tDl 
we identify the origin of the difference on equilibrium en- 
tropy between refs. [2^ andls?! to confirm that the "mem- 
ory corrections" by Ivanov et alr^ are unnecessary. 



II. PERTURBATION EXPANSION WITH 
KELDYSH MATRIX 

A. Contour-ordered Green's function 

We consider an assembly of identical bosons/fermions 
whose total Hamiltonian at time t is given by 

nit)^Ho + H'{t)+Hi,,ta{t). (4) 

Here Hq denotes the kinetic energy, H'{t) is a one-body 
time-dependent perturbation satisfying H'(~oo) — 0, 
-ffint is a two-body interaction, and a{t) is some adia- 
batic factor given by a{t) — 9{—t)e°+* +9{t) , for example, 
with 9 the step function and 0+ an infinitesimal positive 
constant. The system at < = — oo is assumed to be in some 
thermodynamic state described by a density matrix cor- 
responding to Hq. Thus, we here need not consider from 
the beginning the contribution from the path along the 
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imaginary time axis^^^ i.e., "initial correlation," in the 
perturbation expansion with respect to -ffint- 

The explicit expressions of Hq, H'{t) and _ffi„t are given 
in second quantization by 



(5a) 



(5b) 



B. Feynman rules in Keldysh space 

It is desirable to find a simple and compact method 
to carry out the perturbation expansion directly on the 
real-time contour oi —oo<t< oo. This is possible for the 
two-body interaction of eq. ((5c)) (and also for impurity 
potentials), as explained below. We first divide the con- 
tour C into Ci and C2, each running from — 00 to cxd and 
from 00 to —00, respectively. Accordingly, we write the 
integration in eq. ([7]) as a sum of the two contributions: 



Hint = 2 7 '^^^y dW(r-r')V'^(^)V'^(r')V'(^')V'(r) . 

(5c) 

Here m is the particle mass, U is an external potential, 
and V{r) = T^(|r|) is a two-body interaction which can 
be expanded in Fourier series as 



h) 



3 '^1 



(6) 



The spin degrees of freedom will be suppressed for the 
time being. 

We now adopt the interaction representation with re- 
spect to H^°\t) = HQ + H'{t). Then, the time evolution 
of the system is described by the unitary operator: 



Sc = Tc exp 



^ / 77int(i^)dt'^ 



(7) 



Here C is a round-trip contour along the real-time axis 
from t = — 00 towards t = oo^^^ Tq denotes the contour- 
ordering operator along C, and -ffint(i*^) is the interaction 
representation of 7?inta(<) with respect to H^^\t)r-'^-^'^ 

We next introduce the contour-ordered Green's func- 
tion by 



G(l^,2^) 



-(rcV^«(i^)^|,(2^)), 



(8) 



where V''h(1'^) is the Heisenberg operator for ip{ri) with 
1*^ = riti . The perturbation expansion of G with re- 
spect to iJint may be carried out in the same way as 
in the equilibrium theory^^s^li^ by using Feynman dia- 
grams. Indeed, one only needs to change the imaginary- 
time contour of the equilibrium theory into the real-time 
contour C. It hence follows that G satisfies the Dyson 
equation?^ 



G(l^,2^) 

j 3'^)G(3^, 2^) d3^ = 5{l^, 2^) , (9) 



where S denotes the irreducible self-energy. However, 
the round-trip contour C is not convenient for practical 
calculations, since time t appears twice on C with differ- 
ent orders. Thus, we usually have to introduce additional 
contour indices to distinguish them,^^'^^ which make the 
actual calculations rather cumbersome and complicated. 



At- = dt - 

IC J-oo(Ci) 

We next introduce the vector: 
>(li)n 



At. 



= (C2) 



(10) 



^(1) 



V'(i') 



^t(i) ^ [v,t(li) V;t(l2)] , (11) 



where '0(1"') denotes the interaction representation of 
■0(ri) with V = riti'. Then eq. ([7]) can be rewritten 
in terms of '0^, ip and the normal-ordering operator.2£ TV 
as 



Sc = Tc cxp 



-J I'' 



2 ^ 



1') 



x7/;t(i)v^(i)v^t(i')^3^(i') 



(12) 



where 1 = riti with —00 < ti < 00, V is defined by 
V{1 - l') = S{ti - t[)V{ri ~ r[), and is the third Pauh 
matrix. The equivalence of eqs. (O and (|12p may be 
checked easily by writing eq. without using J\f. The 
interaction in eq. (|12p can be expressed diagrammatically 
as Fig. 1 . The expression (|12p is quite useful for our pur- 
pose, because (i) the pairs and ^'^ (l')f3'0(l') 
can be moved around anywhere within the Af and/or Tc 
operators in the perturbation expansion, and (ii) a con- 
traction of with tp^ij) automatically yields a 2x2 
matrix {Tcip{i)'4'\j))^i where the subscript denotes the 
average with respect to H'^^\ Also, the final contrac- 
tion within a closed particle loop can be transformed as 

{Tc^^{i)M{i,j)^{j))o = TVA4(i,j)(rcV^(j)V^t(i))o with 
denoting some matrix product of contractions. 
Now, one may realize that the perturbation expansion 
can be carried out compactly on the real-time axis with- 
out the ambiguity on the limits of time integrations^^ nor 
the complexity from the contour indicesJ^ii^ We intro- 
duce the matrix Green's hmction for this purpose: 



G{1A') ^ -'-{Tci^n{l)A{n 
"Gii(l,l') Gi2(l,l') 

G2l(l,l') G22(l,l') 



(13) 



Here Gi2(l, 1') = T(*A)(0|,(1')0w(1)) = G<(1, 1') and 

G2i(l, V) = -{i/h){,Pn{l)^\i{^'))=G>{l, V) are the cor- 
relation functions introduced by Kadanoff and Bayrtiii 
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FIG. 1: A diagrammatic expression of the interaction in eq. 
(|12|) . The outgoing and incoming lines represent and t/j, 
respectively, whereas the broken line corresponds to the inter- 
action potential. The open and filled circles denote the unit 
matrix and the third Pauli matrix fa, respectively. 



FIG. 2: First-order diagrams for the matrix Green's function. 



with the upper (lower) sign corresponding to bosons 
(fermions). They satisfy 

Gi2(l, 1') = -G^2(l', 1) , G2i(l, 1') - -G;,{1', 1) . 

The diagonal elements can be written explicitly with re- 
spect to the off-diagonal elements as 

Gn(l, 1') - 0(t[-t,)GMh 1') + e(ti-t'i)G2i(l, 1') , 

(14b) 

G22(l, 1') - 0(il-i'i)Gi2(l, 1') + 0(<'i-ti)G2l(l, 1') , 

(14c) 

so that Gii(1,1') = -G;2(1',1) and Gii-fG22 = Gi2+G2_i. 
Thus, there are only two independent components in G, 
i.e., Gi2 and G21. However, all the four elements are 
necessary in the perturbation expansion. Equation (I14p 
can be expressed compactly in terms of G as 

G(l,l') = -fiGt(l',l)fi, TrG-:TrGfi, (15) 



with fi denoting the first Pauli matrix. Equation ([1| 
will be useful later. 

The Fcynman rules to calculate G are summarized as 
follows: (i) Draw all possible nth-order connected dia- 
grams, (ii) With each such diagram, associate a factor 



(»ft)"(±l) 
2"n! 



(16) 



where £ denotes the number of closed loops. Note that 
topologically identical diagrams appear nl times, (iii) 
For each line arriving at 1 from 2, associate the matrix 
G("ni> 2) or f^&^^l, 2), and multiply it from the left of 
the matrix arriving at 2. (iv) If the time arguments of 

G^"^ are equal, we need the replacement Cf^ , G'§2 -^G^y} 
due to the operator M in eq. p2|) . (v) Integrate and 
sum over all the internal variables, and take Tr for every 
closed particle line, (vi) The spin degrees of freedom 
can be included easily by multiplying every closed-loop 
contribution by 25' -1-1, where S denotes the magnitude 
of spin. 

For example. Fig. 2 enumerates topologically distinct 
first-order diagrams for G. The corresponding analytic 



expression is given by 

G(i)(1,1') j A2 j d2'y(2-2'){±G(l,2) 

X [G(2, l')Trf3G(2', 2') + f3G(2, l')TrG(2', 2')] 
+G(1, 2) [G(2, 2')f3 + f3G(2, 2')] G(2', 1')} , (17) 

where we have replaced by G on the right-hand side 
to include the renormalization effects. 

The Dyson equation ([9]) is transformed into a matrix 
form as 



6(1,1') 



dti 2m 

E(l,2)f3G(2,l')d2 = 5(l,l')f3. (18) 



The appearance of f3 between S and G is due to eq. pOj) . 
whereas in front of (5(1, 2) originates from the anti-time 
ordering on G2. Equation (|18p is expressed alternatively 
in an integral form as 

G(l,l') =G(°)(1,1') 

+ / d2 / d2' G(°) (1, 2)f3S(2, 2')f3G(2', l') , (19) 



with 

g(")(i,i') = 



dti ' 2m 



■n ^ 
in— — 



n -1 



-U{1) 



6{l,l')f3. (20) 



Comparing eq. (|17p with the second term on the right- 
hand side of eq. (|19p. we identify the first-order self- 
energy with renormalization as 

SW(1,1') 

- ±'-^S{l, 1') / d2 1/(1-2) [i Trf3G(2, 2) + f3TrG(2, 2)] 



-K^F(l-l')[f3G(l,l') + G(l,l')f3], 



(21a) 



where 1 denotes the 2x2 unit matrix. One can check that 
the off-diagonal elements cancel out in eq. (|21ap . We also 



have to consider the Feynman rule (iv) above. Equation 
(|21ap is thereby simphfied into 

s(i)(l,l')=f3l]H^(l,l'), (21b) 
where is the Hartree-Fock self-energy: 

1') J d21/(l-2)z;iGi2(2,2) 

+V{l-l')ihGi2{l,l'). (22) 

Note E"''(1,1') = [S"''(1',1)]*, i.e., it is Hermitian. 

Thus, the above Feynman rules enable us a straight- 
forward and automatic perturbation expansion of G with 
neither using the contour indices nor worrying about the 
ordering of Tc ■ 

C. ^-derivative approximation 

In carrying out practical calculations, we are almost al- 
ways obliged to introduce some kind of approximations. 
In this context, Bayn3™S, presented an extremely useful 
approximation scheme based on the skeleton expansiouf^ 
i.e., the ^-derivative approximation. The functional 
<& — <&[G] was introduced by Luttinger and Ward as 
part of the exact thermodynamic functional.'^ The $- 
derivative approximation was successively suggested by 
Luttinger— in the equilibrium theory, but has turned out 
to be especially useful for dynamical systemsi^ It has 
the following advantages: (i) it becomes exact if all the 
terms in the skeleton expansion are retained; (ii) various 
conservation laws, which have crucial importance to de- 
scribe dynamical systems, are obeyed automatically; (iii) 
the vertex corrections, or the Landau Fermi liquid correc- 
tions in a different terminology, are naturally included; 
(iv) n-particle (n = 2, 3, • • • ) correlations can be obtained 
with the same approximation scheme, i.e., there is a defi- 
nite prescription here to treat the BBGKY hierarchyji^ A 
detailed study on the dynamical ^-derivative approxima- 
tion has also been performed by Ivanov et aL— for the 
contact interaction. We describe it for the general two- 
body interaction V in terms of the present perturbation 
expansion scheme. It is shown explicitly in AppendijXl 
that various conservation laws are automatically satisfied 
in the ^-derivative approximation. 

Let us define the functional $ in terms of eq. (fT^ by 



<i>EE [(ln5c)o-l] 



skeleton, G<")-^G ' 



(23) 



Thus, $ formally consists of infinite closed skeleton dia- 
grams with replaced by G.— The Feynman rules to 
calculate $ are exactly the same as those of G which are 
given around eq. (|16p . The only care necessary is that we 
have (n — 1)! topologically identical diagrams here. The 
exact irreducible self-energy is obtained from $ by 



E(l,l') = ±f3 



0--0 G> 



FIG. 3: First-order diagrams for $. 



The necessity of 73 on both sides may be realized from 
eq. (fT9|). 

The ^-derivative approximation denotes retaining 
some partial diagrams from the infinite series for $ and 
determining G and E self-consistently by eqs. (fT8|) and 
([24| . It follows from eqs. ^ and ^ that E thus ob- 
tained also satisfies 



E(l,l') = -fiEt(l',l)fi 



(25) 



The first-order diagrams for $ are given in Fig. 3. They 
correspond to 

$(1) J dlj dl'y(l-l') [TrG(l,l)Trf3G(l',l') 

±TrG(l,l')f3G(l',l)] . (26) 

Then E^^' is calculated by eq. ([24]) to yield eq. (|2Ta|) . 

Next, Fig. 4 enumerates topologically distinct second- 
order diagrams. The corresponding analytic expression 
is given by 



$(2) 



222! 



dl / dl' / d2 / d2'y(l-l')T^(2-2') 



X [1VG(1, 2)G(2, l)Trf3G(l', 2')f3G(2', l') 
+TrG(l, 2')f3G(2', l)Trf3G(l', 2)G(2, l') 
±2TrG(l,2')f3G(2',l')f3G(l',2)G(2,l)] . (27) 

The self-energy E*^^) is obtained by eq. ([24]). Expressing 
the result as a single matrix, we observe that the ele- 
ments of E^^-* satisfy exactly the same relations as G^ in 
eq. (HI]), in accordance with eq. ([^5]) . The off-diagonal 
elements are given by (i^j) 

1') = ^{hf /d2 / d2' y(l-2)y(l'-2') 



-T3 . 



(24) 



x[G,,(l,l')G,,(2',2)G,,(2,2') 
±G,,(l,2')Gj,(2',2)Gy(2,l')] . (28) 

The symmetry of S^^^ just mentioned is clearly a general 
property of the higher-order contributions to E, as may 
be checked order by order. Combining it with eq. (|21bp . 
we now realize that En and E22, which satisfy eq. (j25p . 
can be written more specifically as 

En(l, 1') = EHF(i, 1') + 0(t;-ii)Si2(l, 1') 

+0(ii-i;)E2i(l,l'), (29a) 

E22(l, 1') = -E«F(1, 1') + 0(ii-t;)Ei2(l, 1') 

-K0(t'i-<i)E2i(l,l'), (29b) 
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FIG. 4: Topologically distinct second-order diagrams for <1>. 



with S"^ given by eq. 

Besides the one-particle Green's function, the $- 
derivative approximation also provides us with a defi- 
nite and consistent evaluation scheme for higher-order 
correlations, as already shown by Baym and KadanofS^ 
and Baym"^^ on the equilibrium imaginary-time contour. 
However, this topic on the Keldysh contour seems not 
to have been paid due attention in the literature. Espe- 
cially important among them is the two-particle correla- 
tion function: 



/C,,,h(12,34) = (-- ) (Tc^h(iOV'«(3')V'J,(4')V|,(2^')) 



-Gy(l,2)GfeK3,4) 



(30) 



Note that the arrangement of the space-time arguments 
in K. is different from that of Baym and Kadanoff;^^ the 
present one may be more convenient when regarding K. 
as a matrix. 

To consider higher-order correlations in a unified way, 
we introduce an additional perturbation on the Keldysh 
contour caused by the non-local one-body potential 
W{1^ ,2'-^). Adopting the interaction representation in 
terms of H in eq. ([4]), the time evolution due to W is 
described by the operator: 



S' = Tc exp 



dl / dl>|,(l')f3W^(l',l)f3V'^^(l) 



(31) 

where fa's originate from the transformation pop . The 
Green's function is now given by 



G(l, 1'; W) = ~-{TcS'M^)Ai^'))l{TcS') . (32) 



The correlation function (j30p is then obtained by 



my,{2\2) 



(33) 



where ± is due to the commutation relation between ■0^ 
and "0, and (— l)-'"'"-' cancels the contribution of fs's in eq. 



((3T|) . To calculate K self-consistently in the ^-derivative 
approximation, we start from the Dyson equation (jlSp . 
It reads symbolically as G^^G = 1, where G"^ is given in 
the present case by G~^—f'i[ihdt — K — \J) — f-j,{yV^Y^f'i 
with K denoting the kinetic-energy operator. Therefore, 
the first-order change 8G obeys 5G = G[—&G~^)G . Let 
us substitute -6G-^ =fz\m + 8G){bG I m\m\f-i 
into the first order equation and divide it by SW . We 
thereby obtain an integral equation for eq. (|33p as 

K^^jy (11', 22') = ±Gyv (1, 2!)G,v (2, l') 

^i^Y^ J'isJ'i'i y"d3'y"d4'G.fc.(l,3')Gfc..(3,l') 

kk'U' 

xTkk'.w (33', 44')/C„.,jy (44', 22') , (34) 
where F is the irreducible vertex defined by 



r..,,,,,(ll',22')^T^(-l^^+'' 



' 6Gjj> {2,2') 
i 

"ft5G,,-(l,l')<5Gjy(2,2') 



(35) 



We have used eq. (124)) to derive the second expression 
of eq. pS)) . Thus, once $ is given explicitly as a func- 
tional of G, the two-particle correlation ([50]) can also be 
calculated by eqs. ([M)) and ([55)1 . 

The integral equation (p4|) may be solved iteratively to 
obtain a formal solution: 

^ = ±(iT*ftGGf)"^GG, (36) 
where 1 and GG are matrices defined by 

(i),,.jy(ll',22') = %5,,,.5(1,2)<5(1',2'), (37a) 

(GG),,^,y(ll',22') -G,,v(l,2')G,,,(2,l'), (37b) 

respectively. 

Equation (|35p clearly has the symmetry: 



r,,,,,y(ii',22') -rjy,,,,(22',ii'). 

It also follows from eq. (fTSl) that Tar jj, satisfies 
Tu'jj'in',22') = - ifi)ik{fi)k'i'{fi)jiifi)i,j 



kk'U' 



[Ffe,,,n(l'l,2'2)]* 



(38b) 



The expressions of f|^-*^^-, (11', 22') for ti = 1, 2 are given 
in Appendij^B] to see their structures explicitly. 



D. Keldysh transformation 



As seen from eq. ([T4| . the four elements of G are not 
independent. This redundancy in G is removed by the 
following modified Keldysh transformation] 



G^ = Lt-^GD ^ 



G^ GK 
G^ 



(39) 
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where L is defined by L = -^(i — 1x2) with f2 denoting 

the second Pauli matrix. Thus, the 21 element of 
vanishes, and the others also satisfy 

[G"(l, 1')]* - G^(l', 1) , [G^(l, 1')]* = -G^(l', 1) . 

(40a) 

We hence realize that G^ and G^ form an alternative set 
of two independent elements in the Keldysh space. In- 
deed, they are connected with the Kadanoff-Baym func- 
tions G12 and G21 as 

G«(l, 1') = 0(ii-i'i)[G2i(l, 1')-Gi2(l, 1')] , (40b) 



GK(1,1') = Gi2(1,1') + G2i(1,1') 



(40c) 



It also follows from cqs. ([25]) and ^ that t^ = Lf:itD 
can be written as 







(41) 



Its elements satisfy 

[E«(1,1')]*=S^(1',1): 



P^(l,l' 



= -s^(iM). 

(42a) 

The quantities E^ and are given in terms of E^^, 
E12 and E21 as 

E«(l, l')-S"^(l, l') + ^^(ti-i'i)[S2i(l, 1')-Si2(l, 1')] , 

(42b) 



eK(1,1') = Sl2(l,l') + S2l(l,l')• 



(42c) 



It is worth pointing out that the Keldysh transformation 
is not useful in the perturbation expansions of G and 
$, since it obscures the basic symmetry of them. The 
transformation should be carried out only after finishing 
the expansion in terms of G. 

Applying the Keldysh transformation to eq. p8)) . we 
obtain the Dyson equation for G^ as 



GK(1,1') 



- / EK(l,2)G'^(2,l')d2 = (5(l,l')i 



(43) 

Thus, the equation for the retarded function G^ is com- 
pletely decoupled from that of the Keldysh component 
G^. 

Finally, it follows from eqs. jHSl), ([SHI) and (|41|) that 
the variation (5E^ is connected with 5G^ as 

6Y}^,{l'A)^±zhY, /d2| d2'r,^,_^.^.,(ll',22') 

jj' 

x<5Gf,,(2,2'), (44) 



with r^, ■■, defined by 



r^,_^.^.,(ll',22') ^ {Lf^U{L^),.ALf,),i{P) 



kk'W 



crfcfc,,„,(ll',22'). 



(45) 



We realize from eq. ((38|) that F^, jj, satisfies 

fK _^.^.,(11',22') = Ff^.,,,,,(22', 11') , (46a) 



FK ^^.^.,(11',22') = - ^ {if2Mtf2)kH'{tf2)Air2)l'y 
kk'W 

x[FK,,^(1'1,2'2)]*. (46b) 

The quantity the advantage that we only need 

to consider its 9 elements instead of 16 in Taijji due to 
the vanishing 21 elements in E^ and G^. 



III. GRADIENT EXPANSION 

The theoretical framework of ijll] enables us a formally 
exact microscopic treatment of nonequilibrium dynami- 
cal systems. However, the coupled equations ([24| and 
(|43| are still too difficult to solve practically. It is desir- 
able to reduce their complexity down to a tractable level 
without loosing the physical essentials. To this end, we 
here adopt the Wigner representation and subsequently 
carry out the gradient expansion to eqs. and (|43p . 
To be specific, the Wigner representation of G(l,2) is 
defined through 



G(l,2) 



^^G(p£, rl2^l2)e^(^•'=--^*-)/^ (47) 



with_ri2 = (ri+r2)/2, ti2 = {ti+t2)/2, ri2 = ri-r2, 
and ti2 = ti—t2- The gradient expansion denotes an ex- 
pansion with respect to rt of G{pe, rt). It forms a well- 
established basis for the microscopic derivation of various 
transport equations such as the quantum^'' and classical 
Boltzmann equation. The basic assumption is that the 
scales of the space-time inhomogeneity are much longer 
than the microscopic scales to achieve the local equilib- 
rium such as the mean-free path and the collision time. 
This condition is well satisfied in most of nonequilibrium 
steady states without time evolution such as Rayleigh- 
Bcnard convection4^^>27^iS9^ 



A. Spectral and distribution functions 

As seen in eq. (HH), there are essentially two indepen- 
dent components in G, i.e., G12 and G21. We here in- 
troduce an alternative pair of independent components, 
i.e., the spectral function A and the distribution func- 
tion 0,^^'^^ which turn out to be more convenient. The 
spectral function A is defined by 



^(l,2)EEj[G2i(l,2)-Gi2(l,2)] 

= ^(V'w(1)^|,(2)t^|,(2)Vh(1)). (48) 

Let us expand A as eq. (|T7)) . It then follows from the 
equal-time commutation relation of the field operators 
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that A{pe, rt) satisfies the sum rule: 

/oo 
A{pe, rt) de = 2n. 
-OO 



(49) 



We also conclude from A*{1,2) = A{2, 1) that A{pe,rt) 
is real. We next introduce the distribution function 
directly in the Wigner representation as 

(j){pe,ri2ti2) 

/I 
d^ria - (V-l, (2)^„ (l))e-*(P-^--^*-)/'"' 



A{pe,ri2ti2) 



(50) 



We find from A* (pe , rt) ^ A{pe , rt) and (?/'|^(2)?/'-h(1)}* 
= (?/'|^(1)V'-h(2)) that 0(p£, rt) is also real. Thus, both A 
and are real in the Wigner representation. They form 
an alternative pair of independent quantities in G. 

In the equilibrium theory, (/) is just the Bose/Fermi 
distribution function so that we only need to calculate 
A. For nonequilibrium systems, in contrast, we have to 
determine A and </> simultaneously. The derivation of 
the quantum transport equatioi*i^ii^ii^iiiii^ii^i2ai2ii^i2^ 
amounts to integrating out the spectral function A, which 
contains detailed information on the density of states, to 
obtain a single equation for 0. However, we shall proceed 
without the approximation here. 

B. Wigner representation of G and E 

The Wigner representation of G can be written in 
terms of A and </). Indeed, the off-diagonal elements are 
transformed into 



Gi2{pe,rt) = ^iA{p£,rt)4>{p£,rt) 



(51a) 



G2i{pe,rt) = ~iA{p£,rt)[l±(i){pe,rt)]. (51b) 

Using eqs. (PH)) and ([5T|) . we also obtain the Wigner rep- 
resentations of G^ , G^ and G^ as 



G''ipe,rt)^[G^{pe,rt)r= J — 



de' A{pe\rt) 



e+ ~ e' 

(52a) 



G^{pe,rt) = ~iA{pe,rt)[\ ±24>{pe,rt)] , (52b) 

with £_|_ =e-l-iO+. Since it may not cause any confusion, 
we use the same symbols Gij in both the coordinate and 
the Wigner representations. All the symbols given below 
without arguments belong to the Wigner representation. 

Let us move on to the self-energies of eqs. ((24)) and (|4T|). 
Following Ivanov et al.^ we express the Wigner repre- 
sentations of their independent components in a form 
similar to the above expressions as 



T.i2(.P£,rt) = ^iAY,(pe,rt)(i)Y,(pe,rt) 



(53a) 



S21 (pe, rt) = As (pe, rt) [1 ± {pe, rt)] , (53b) 



and 



=sHP(p.«,+ r±:i£(E^. ,54.) 



S^(pe,rt) -iA^(pe,rt)[l ±2(t)^{pe,rt)] . (54b) 



It follows from eq. ([25)1 that As and 0s are also real. 
They are functionals of A and in the ^-derivative ap- 
proximation. 

The operator in the square bracket of eq. ([15)) can also 
be transformed into a matrix of the space-time coordi- 
nates by multiplying it by 5(1,3) from the right. Its 
Wigner representation is given by 



Gq ^iP£,rt) = e 



V 

2m 



U{rt) . 



(55) 



C. Gradient expansion 



We now consider the matrix product of the space-time 
coordinate. It is transformed into the Wigner represen- 
tation as- 

dh2 /"d^fia e-^(P-^i^-^*i^)/'^ f d3C(l,3)Z?(3,2) 

-00 J J 

^ C{pe,ri2ti2) D{pe,ri2ti2) , (56) 
where the operator (g) is defined by 
C{pe,rt)(E)D{pe,rt) 

ih 



= exp 



[dr -dp' - dtde' - dp dr' + dedf 



xC{pe,rt)D{p'e',r't" 



(57a) 



p' —p.e' —e ,r' —r ,t' —t 



with dr = d/dr and dt = d/dt. The identity (|56p can 
be proved as follows: write G and D on the left-hand 
side in the Wigner representation; expand ri^ti^ and 
''32^32 from ri2ii2; remove risiis and f32f32 by using 
f^^ggV '^ia/?! = ^ifi^Qip' ■ f 13 /n.^ Q^^ . perform partial 

integrations over internal momentum-energy variables; 
carry out the integrations over 3 and ri2ii2. The first- 
order approximation to eq. ()57ap yields 



C(^DkGD+'-^{C,D}, 



(57b) 



where the curly bracket denotes the generalized Poisson 
bracket: 



{G,D} 



dC dD dC dD dG dD dG dD 
dr dp dt de dp dr de dt 



(58) 



Finally, it should be pointed out that both the Wigner 
transformation ()47p and the gradient expansion (|56p need 
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essential modifications in the presence of the electromag- 
netic field. Here we need a special care on the gauge 
invariance of the equations in order to appropriately ob- 
tain (i) the Hall terms2ii22i2^ and (ii) the pair poten- 
tial of superconductivity as an effective wave function of 
charge 2e.— To be specific, the Wigner transformation 
(ITf)) should be defined as eq. (7) of ref. [2l| or eq. (21) of 
ref. [23l. an d eq. above has to be replaced by eq. (36) 
of reflH. 



D. Dyson equation in the Wigner representation 

We now transform the Dyson equation into the Wigner 
representation within the first-order gradient expansion. 
Following Keldysh/^ we start from eq. (|43)) rather than 
eq. (fTS]) adopted by Kadanoff and Baymii and Ivanov 
et alM' Indeed, eq. has the advantages that (i) its 
21 element vanishes and (ii) the 22 element is complex- 
conjugate of the 11 element. Hence we essentially need to 
consider only the first row of eq. (|43|) . which completely 
determines the two independent components of G, i.e., A 
and (j). Thus, we can see the structure of the equations 
more clearly in the present approach. 

Using eq. dSH) and (|57b|) . the 11 element of eq. ^ is 
transformed into 

(Go-i-I]«)G« + |{Go-i-SR, G«} = 1 , (59a) 

where Gq^ is defined by eq. ([55]) . The replacement R— * A 
in the superscript yields the equation for the 22 element. 
Taking its complex conjugate and noting eq. (|52ap . we 
have an alternative equation for G^ as 

(Goi-S^)Gi^- |{Goi-S«,Gi^} = 1. (59b) 



Let us add eqs. (|59ap and (j59bp . Wc then obtain 



= (Go 1 - S^^ ) 



Ra-1 



(60) 



We realize from eq. (|52ap that this is the equation to de- 
termine the spectral function A for a given S^. Since the 
11 and 22 elements of eq. (|43p were equivalent before the 
gradient expansion, it is natural to ask whether this prop- 
erty is still retained between eqs. (|59ap and (|59bp . To an- 
swer the question, let us subtract eq. (|59b|) from eq. (|59ap 
and substitute eq. (|60p in the resulting equation. We then 
obtain = {{G^)-\G^} = ~{G^)-^{G^,G^}, which 
is just = 0. We have thereby confirmed the equivalence 
between eqs. (|59al) and (j59bp . Wc now realize that A 
can be determined locally by eq. (j60p without space-time 
derivatives within the first-order gradient expansion. 

It follows from the retarded nature of G^(l,2) in eq. 
(I40bp that all the singularities of G^ (pe, rt) in eq. ((SD)) 
lie on the lower half of the complex e plane. This implies 
ImE^ (pe, rt) < and hence ImG^(p£, rt) < 0. Using eq. 
(|52ap . the latter condition can be written alternatively 
as 



Next, the gradient expansion to the 12 element of eq. 
(gSl) leads to 



(Go 1-SI^)GK - S^^G^ + ^{Go-^-S«, G^} 
-|{EK,Gn = 0. 



(62a) 



Taking its complex conjugate and using eqs. (|52p and 
([5l)) . we have an alternative expression of eq. (|62ap as 



(Go-I-S^)GK - S^G^ - ^{Gq-I-S^, G^} 
+ |{E^G«} = 0. 



(62b) 



Let us subtract eq. (|62bp from eq. (|62ap . substitute eqs. 
([52ll and ([54]) into the resulting equation, and use {Gq^— 
S" , G^ } = 0. We thereby obtain 



AAy 



{Go'-ReS^,^0} - {As0E,ReG«} 



The left-hand side of the equation consists of terms with 
space-time derivatives, whereas the right-hand side de- 
notes the collision integral which vanishes in equilibrium. 
Hence it is appropriate in the first-order gradient expan- 
sion to replace ips on the left-hand side by (j). The ap- 
proximation was originally suggested by Botermans and 
Malfiiet^° and adopted explicitly by Ivanov et al}^ The 
idea also has a close relationship with the Enskog series 
for solving the Boltzmann equation , ^'^i-'^^ i.e., the expan- 
sion from the local equilibrium. The procedure leads to 

{Go"^-ReI]^,^0}-{As0,ReG^} = C, (63) 

where C denotes the collision integral defined by 

AAs(0s - (j) 



C = 



G21S12 — G12S2I 

T ; • (64) 



Equation (|63p determines cj) for given A and E. 

Following Ivanov et al.^ we now ask the question of 
whether eqs. (j62ap and (|62bp still hold the equivalence 
which was present before the gradient expansion between 
the 12 element of eq. p3)) and its complex conjugate. Let 
us add eqs. (|62ap and (j62bp . After the same procedures 
as above for the subtraction, we obtain 

{A^(t>,A}^{A^,A(l>} 
4 

_ As<?!)sReG^ - (Go ^-RcER)A0 



A(pe,rt) > 0. 



(61) 



This equation is identical with eq. (j63p . Indeed, multi- 
plying it by Ae/(Go"^— ReS^) yields eq. ([63|l in disguise. 
This may be seen more explicitly by (i) expressing ReG^ 
and A in the two apparently different equations with re- 
spect to ReG^~^ =Go ^— ReE^ and A-^, and (ii) trans- 
forming the gradient terms into derivatives of ReG^~^, 
Ay: and (pi^^ Thus, with the approximation (p-^—Kp in the 
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gradient terms, eqs. (|62ap and (j62bp recover the original 
equivalence. 

Equations ([60|) and ([63|) form coupled equations to 
completely determine the real quantities A and (f> for a 
given E. Indeed, eq. (j63p is real, and the real and imag- 
inary parts of eq. ([50]) are connected by the Kramers- 
Kronig relation. 

In AppendijCl we derive basic conservation laws in the 
first-order gradient expansion of the ^-derivative approx- 
imation. 



E. Self-energy in the Wigner representation 

The original ^-derivative approximation denotes solv- 
ing eqs. ((24|) and (j43|) self-consistently. Having performed 
the first-order gradient expansion to the Dyson equation 
(I43p . we also have to specify a consistent approximation 
scheme to the other equation ((24|l . This issue seems not 
to have been given an explicit consideration before, how- 
ever. 

As noted below it, eq. ([60)1 is correct up to first order 
in the gradient expansion. This implies that the local 
approximation to S'^ is sufficient for solving eq. (|60p . As 
for eq. ([55]) . all the terms on the left-hand side include 
space-time derivatives, and the collision term of the right- 
hand side is connected by equality with the left-hand 
side. Thus, eq. (|63p is a first-order equation where every 
quantity should be evaluated locally without derivatives. 
With these considerations, we now conclude that we have 
to apply the local approximation to eq. ([M)) . We have 
thereby reached the definite prescription to evaluate S in 
terms of G, so that eqs. (|60l) and ((63)l now form closed 
equations for A and (j). 

Using eq. ([5]), the Hartree-Fock self-energy ([^ is 
transformed into the Wigner representation as 



j^^iVo ± W)Gi2(p'£', rt) . 



(65) 



Also, the local approximation to eq. ([28l) yields [ij^j) 



n'^ d^PkdSk 1, 



k=2 



X {2TThY5 {p+P2 -P3 -P4)^ (£ + £2 -£3 -£4) 

xGji{p2S2,rt)Gij{p3e3,rt)Gij{p4e4,rt) . 



(66) 



Writing G12 = TiG<, G21 = -iG> , S12 = TiS< and 
S21 = — in the Wigner representation, we find that 
eq. (|66)) is identical with eq. (4-16) of Kadanoff and 
Baymj^ as they should. 

Contrary to the local approximation adopted here, 
Ivanov et alr^ emphasized the importance of consider- 
ing the first-order gradient corrections to the collision 
integral, which they call memory corrections. Their mo- 
tivation towards this conclusion seems to be stemming 
from the expression of equilibrium entropy obtained by 



Carneiro and Pethickj^ see the paragraph at the end in 
§5.4 of Ivanov et alJ^ Indeed, they have shown that the 
Carneiro-Pethick expression cannot be reproduced from 
their dynamical equations without the memory correc- 
tions. On the other hand, we have already seen that the 
local approximation should be sufficient for the collision 
integral. Thus, one may wonder which statement is cor- 
rect and where the discrepancy originates from. We will 
show that: (i) the Carneiro-Pethick expression of equi- 
librium entropy is not correct due to an inappropriate 
treatment of energy denominators in their calculation of 
the thermodynamic potential; and (ii) the local approx- 
imation without the memory corrections leads to an ex- 
pression of dynamical entropy which is compatible with 
the equilibrium expression. Thus, the memory correc- 
tions should not be incorporated within the first-order 
gradient expansion. 



IV. ENTROPY 

We are now ready to discuss entropy in noncquilibrium 
statistical mechanics. We first derive an equation of mo- 
tion for entropy density in ijIV Al We then prove the 
i?-theorem for a limiting case in i^IVBl Finally, we pro- 
pose a principle of maximum entropy for noncquilibrium 
steady states in i^IV CI 



A. Equation of motion for entropy density 

Let us multiply eq. (|63p by fee ln[(lzb^)/(/)] and integrate 
it over p and e. We next write fce ln[(l ± (p) / (l)]d(j) = da 
in the resulting equation, where a denotes entropy of the 
noninteracting system)^ 



o- = fcB[-01n(/) ± (1 ± (f)) ln(l ± <j>)] 



(67) 



Thus, the derivatives of (p on the left-hand side are trans- 
formed into those of a. We then perform partial integra- 
tions over p and e. We thereby arrive at an equation of 
motion: 



d^p de ^ , 1 ± I 
C in — — 



(27rft)4 

with s = s{rt) and js=3s{i't) defined by 



(68) 



= h 



d^p de 
(27rft)4 

d^p de 



A 



9(^0 ^-ReS^) 



■A, 



dRcG^ 
de 



(69) 



-A 



9(^(7^ -RcS^ 
dp 



A- 



aReG« 



dp 



(70) 

Let us explain the quantities in these expressions once 
again for an easy reference. The quantity C is the colli- 
sion integral (j64p . cf) the distribution function defined by 
eq. ([5D|) . A the Wigner transformation of eq. called 
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the spectral function, Gq ^ defined by eq. ([55)) . and the 
retarded functions and given by eqs. (|52ap and 
()54ap . respectively. The basic quantities A and (j) should 
be determined self-consistently by eqs. ([SO]) and us- 
ing the self-energy of the local approximation, as dis- 
cussed in detail in §111 El 



The right-hand side of eq. ((68)) denotes net change of 
entropy due to collisions, whereas s and js on the left- 
hand side can be regarded as the entropy density and 
the entropy flux density, respectively. Indeed, eq. ([M]) 
agrees completely in form with the equilibrium expression 
of entropy, i.e., eq. (5) of ref. [2^4^ This may be shown 
explicitly from the latter by: (i) noting d(j>/dT — —da/de 
for (j)— [e'^^''^'^ ±1)^^ in equilibrium; and (ii) performing 
a partial integration with respect to s. Thus, eq. (j69p 
is compatible with equilibrium statistical mechanics and 
may be regarded as an expression of the nonequilibrium 
entropy density. 



Equation ([55]) was obtained by Ivanov et al.^^ as the 
equation for "Markovian entropy flow." However, they 
claim that there is additional contribution to entropy 
called "memory effects," which originates from the gradi- 
ent terms in the collision integral on the right-hand side 
of eq. ((63)) . Note that the gradient terms in the collision 
integral have been discarded in the present formulation 
with rationales given in the second paragraph of ijlllEI 
By including the memory effects, Ivanov et ali^ could 
obtain an expression of entropy which is compatible with 
the equilibrium entropy derived earlier by Carneiro and 
Pethick,^ 

However, Carneiro and Pcthick'^^ obtained their ex- 
pression with the zero-temperature time-ordered Gold- 
stone technique where there may be ambiguity as to how 
to deal with vanishing energy denominators. Indeed, 
Carneiro and Pethick added an infinitesimal imaginary 
quantity to every energy denominator in their calculation 
of $ in equilibrium. They thereby found a contribution 
to $ from singularities in the energy denominators, i.e., 
the "on-energy-shell" term; see the arguments in §IV.A 
of their paper It is this on-energy-shell contribution to 
$ which brings the difference between eq. (|69p and the 
Carneiro-Pcthick expression. 

Since the issue has a crucial importance to the whole 
theory, we have reexamined in Appendi jP] whether such 
singular contribution to $ is really present or not. To this 
end, we adopt the finite-temperature formalism of using 
the Matsubara Green's function. A definite advantage of 
the present approach is that no additional regularization 
procedure is necessary. It is thereby shown that the on- 
energy-shell contribution is absent. Thus, it is eq. ([69]) 
which is compatible with the equilibrium expression of 
entropy. The conclusion also implies that we need not 
consider the "memory effects" of Ivanov et ali^ 



B. Entropy production and the //-theorem 

As it has already been mentioned, the right-hand side 
of eq. (l68|) expresses net change of entropy due to colli- 
sions with C given by eq. (|M|) . We hence put 



d'^p de 1 zh I 

T kB I ^-.^ (G21S12 - G12S21) In — — 



(27rft)4 



^5col 

dt 



(71) 

and study this term more closely. It is shown shortly 
below that eq. (|7ip is positive within the second-order 
perturbation expansion. 

Let us substitute eq. ((66)) into eq. ((TT)) and rewrite it 
in terms of the Kadanoff-Baym functions defined in the 
Wigner representation by^^ G^ = zLiGi2 = A(j> and G-* = 
iG2i = A(l± 0). Equation ((TTj) is thereby transformed 
into 



d. 



.(2) 



dt 



X {2TTh)^S{pi +P2 -P3 -P4)S{ei+e2 - £3 - £4) 

/^</^> /^</^> 

[r<<r'>r<<r<> r<> n<n> r<<\ i„ 1 2 ^3 '-^4 
^ L*-^! '^2 ^3 "-^4 ~ '-'I '^2 "-^3 ^-^4 J n> n< n> n< ' 

•-"I ^2 ^3 "-^4 

(72) 

with G^ = G^{pjej,rt) and G^ = {pjSj.rt). Using 
the inequality {x — y)\Ti{xly) > which holds for any 
positive X and we then conclude 



ds'f^^{rt)/dt>{). 



(73) 



Thus, entropy increases by collision in the second-order 
perturbation. This is quite a strong statement in that 
the inequality holds even locally. 

The iJ-theorem is relevant to the space integral of eq. 
over the whole systemi^ii^ 

j s{rt)d\ + j js{rt)-dS^^^ j Scoi(rOdV, (74) 



d 
'dt 



where dS denotes the infinitesimal surface element. We 
review it here for a later extension: Consider an isolated 
system where the second term on the left-hand side van- 
ishes. If the time evolution of entropy by collision is 
globally positive, i.e.. 



^ j scoi(rt) dV > 0, 



(75) 



we obtain the law of increase of entropy for the relevant 
system as 



— I s{rt)d^r > 0. 
dt J 



(76) 



It hence follows that entropy takes its maximum value in 
equilibrium of an isolated system. 

There seems to have been yet no explicit proof for 
eq. ()75p beyond the second-order perturbation. Indeed, 
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the iJ-theorem has been discussed almost exclusively in 
terms of the Boltzmann equation with the two-particle 
(i.e., second-order) collision integral. A general expres- 
sion corresponding to eq. (|72p has been provided for the 
contact interaction as eq. (5.11) of Ivanov et al.^^ A suf- 
ficient condition for the general if-theorem to hold is 
Rm.m' > in their expression at every order of the per- 
turbation expansion, which seems nontrivial to prove at 
present and we defer it for a future study. However, it 
should be noted that Nature clearly adopts the inequality 
()75p beyond the second-order, since no explicit violation 
of the second law of thermodynamics has been reported 
even in systems with strong correlations. See also the 
numerical work by Orban and Bellemans on the time 
evolution of entropy of an isolated system4i 



C. Maximum entropy in nonequilibrium steady 
states 

We now extend the above principle of maximum en- 
tropy to nonequilibrium steady states without time evo- 
lution. Wc consider specifically those cases where there 
is influx of current Jz and/or energy current J^z through 
one boundary perpendicular to the z axis. It follows from 
the conservation laws that there is the same amount of 
currents flowing out through another boundary in the 
steady state; hence these quantities can surely be adopted 
as additional variables of entropy to specify the system. 

First of all, we assume that eq. (TTS]) also holds in steady 
states. This can be proved explicitly within the second- 
order perturbation as eq. ((73)) . Next, we note that there 
is no net inflow or outflow of entropy in steady states 
as well, i.e., the second term on the left-hand side of eq. 
([7^ vanishes. Hence it is quite reasonable to expect that 
eq. ([76]) also holds for steady states under appropriate 
conditions. The question then arises: what are the pa- 
rameters that have to be fixed? In this context, we note 
that the principle of maximum entropy in equilibrium can 
be stated without the magic word of "isolated system" as 
Scq{E, V, N) = maximum, where energy E, volume V and 
number TV are all mechanical variables. This is quite nat- 
ural, since "probability" can only be introduced at first in 
terms of mechanical variables. In contrast, temperature 
T, pressure p and chemical potential fi are all equilib- 
rium thermodynamic variables defined in terms of Scq 
by partial differentiations, i.e., there are no definite defi- 
nitions for them in nonequilibrium. Hence the latter can- 
not specify nonequilibrium states of the system. These 
considerations indicate that we should choose mechani- 
cal variables as independent variables of nonequilibrium 
entropy. We hence add Jz and J^z, which are conserved 
and can be calculated mechanically, as independent vari- 
ables of entropy in the present context. Now, we extend 
the principle of maximum entropy as follows: 

Principle of maximum entropy for steady states: The 
state which is realized most probably among possible 
steady states without time evolution is the one that 



makes S{E, V, N, Jz, J^z) maximum. 

The validity of the principle can only be checked by its 
consistency with experiments. In the next paper we shall 
test it on Rayleigh-Benard convection of a dilute classical 
gas which is typical of nonequilibrium steady states with 
pattern formation4*2^i2&i21i2ai2ai22 It will be shown that 
the convection indeed gives rise to an increase of entropy 
over the value of the heat conducting state. 

Finally, it may be worth pointing out that, once 
S{E,V, N, Jz, Jez) is given explicitly, we may perform 
successive Legendre transformations to change indepen- 
dent variables, just as in the equilibrium theory. 



V. SUMMARY 

We have performed a theoretical study on entropy in 
nonequilibrium statistical mechanics by specifically con- 
sidering an assembly of identical bosons/fermions inter- 
acting via a two-body potential. First, we have presented 
an expression of nonequilibrium entropy density as eq. 
(f69l) . which obeys the equation of motion ([68|l . Thus, 
we can now trace time evolution of entropy in the many- 
body system. Second, we have proposed a principle of 
maximum entropy in mV CI for nonequilibrium steady 
states without time evolution. The validity of the prin- 
ciple will be checked in the next paper by calculating 
the entropy change of a dilute classical gas through the 
Rayleigh-Benard convective transition. 

A conventional theoretical starting point to nonequilib- 
rium systems has been some phenonicnological determin- 
istic equations connected closely with the conservation 
lawsi^ One then performs the linear stability analysis and 
derives some effective equations near the instability point 
such as "amplitude equations" or "phase equations." In 
some fortunate cases one may further be able to construct 
a Lyapunov function from those differential equations;^ 
Note however that this approach is essentially of mechan- 
ical character, as it is completely irrelevant to the con- 
cept of probability. In contrast, little attention seems 
to have been paid to entropy in nonequilibrium systems, 
even near instability points, due partly to the absence of 
an explicit expression of nonequilibrium entropy. Since 
entropy is the key concept of equilibrium thermodynam- 
ics and statistical mechanics embodying "probability," it 
will be well worth studying entropy of nonequilibrium 
systems and their "phase transitions," which will shed 
new light on the phenomena. The theoretical framework 
proposed here may provide a starting point for those in- 
vestigations. Its obvious advantage over the approach 
of the nonequilibrium statistical operator by Zubare"\«^ 
is that one can treat nonequilibrium systems which are 
globally far away from equilibrium, though not locally. 

Homogeneity /additivity has played a key role in con- 
structing equilibrium statistical mechanics. In contrast, 
the present approach may be regarded as an attempt 
to treat open inhomogeneous systems by an extremum 
principle with considering the boundary conditions ex- 
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plicitly. Once it is established that steady states are 
identified correctly with the extremum principle, it will 
be a straightforward task to develop the linear-response 
theory around it in the same way as in the equilibrium 
theoryi^ 
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APPENDIX A: CONSERVATION LAWS 

We here show in terms of the present real-time per- 
turbation expansion in the Keldysh space that various 
conservation laws are automatically satisfied in the 
derivative approximation. We follow essentially the pro- 
cedures by Baym^^ with a slight modification appropriate 
for the real-time contour. 



1. Identities 

To start with, we derive several identities which form 
the basis for proving the conservation laws. 

First, consider the following gauge transformation: 



G(2,l) 



G(2,l)e- 



with 



x(i) 







(Ala) 
(Alb) 



It yields a first-order change in G as 

66(2, 1) = I [x(2)G(2, 1) - G(2, l)x(l)] . (A2) 

However, $ is clearly invariant through eq. (|Aip . With 
eq. (|24p . this invariance of $ reads 



dl / d2 Trf3E(l, 2)f3(5G(2, 1) = . 



(A3) 



Substituting eq. (|A2[) into eq. (jA3p and using is 
arbitrary, we obtain 



/ 



d2 Tr ^ [f3E(l, 2)f3G(2, 1) - G(l, 2)f3S(2, l)f3] 
= . (A4a) 
We can further transform eq. (|A4a[) into 

J d2 [S«(l, 2)Gi2(2, 1) + Si2(l, 2)G^(2, 1) 

-GR(1, 2)S]i2(2, 1) - Gi2(l, 2)2^(2, 1)] = , (A4b) 



where we have used eqs. (jl4bl) . (I29ap . (j40bp . (j42bp and 
the Feynman rule (iv) around eq. (jl6p . This is the basic 
identity obtained with respect to the gauge transforma- 
tion. 

Second, consider the following Galilean transforma- 
tion: 



G(2,l) 
with 



exp 



G(2,l)exp[Vi-fi(ii)J , 
(A5a) 



R{t) 




(A5b) 



It yields the following first-order change: 

(5G(2,1) ^i?(i2)-V2G(2,l) + %i)-ViG(2,l). (A6) 

However, $ is invariant through eq. (|A5P so that eq. (|A3p 
holds also in this case. Substituting (jA6P into eq. (|A3P 
and using R(t) is arbitrary, we obtain 



yQ(i)dVi =0, 



(A7) 



with Q(l) defined by 



Q(l) = Titi- 



d2 [si^(l,2)Gi2(2,l') 



+Ei2(l, 2)G^(2, 1') - G^'il, 2)Si2(2, l') 
-Gi2(l,2)E^(2,l')],,^,. 



(A8) 



Here terms with the self-energy derivative are due to par- 
tial integrations. Equation (|A7p is the identity obtained 
in terms of the Galilean transformation. 

We finally consider the change of time on Gi: t ^ 
6{t)=t + (p{t). Accordingly, G is transformed as 



G(2,l) -^U{t2)Gir292,riei)Uih), 



with 



(d6l/dt)i/4 




(A9a) 



(A9b) 



The factor {d9/dty^^ cancels the Jacobian for t — > 6* in 
eq. p2p . thereby keeping $ invariant in form. We hence 
conclude that eq. (|A3p also holds in this case. The first- 
order change in G is given explicitly by 



,SG,,(2,l) = |<5,i 



4 

'^'(tl) 



Lp{tl 



dti 



Gj,(2,l). (AlO) 



Substituting eq. (jAlOp into eq. (|A3p and using ip{t) is 
arbitrary, we obtain 



d(gint(^l)) 

dti 



- /g,(i)dVi 



(All) 
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where (-ffi„t(ii)) and are defined by 

(i?i„t(ti)) = ±j /d'^i/d2 [S^(l,2)Gi2(2,l) 

+Ei2(l, 2)GA(2, 1) + G^(l, 2)Ei2(2, 1) 
+Gi2(l,2)5]A(2,l)] , (A12) 



ge(l) = iT^h)^J d2 [E«(l', 2)Gi2(2, 1) 

+Si2(l', 2)GA(2, 1) + G«(l, 2)Ei2(2, 1') 
+Gi2(l,2)SA(2,l')],,^,, (A13) 



respectively. Equation (jAlip is the identity obtained 
from eq. 

It should be noted that (i?int(ii)) defined above is ex- 
actly the interaction energy of the system. This can be 
seen as follows: With eqs. (l40b| and (|42b)) . 

the 12 element of eq. (fT8|) is transformed into 



dti 2m 
Ei2(l,3)G'^(3,2)] d3 = 



Gi2(l,2)-y[Ei^(l,3)Gi2(3,2) 

(A14) 



We may alternatively derive an equation for Gi2(l, 2) — 
=F(z//i)(^/'|^(2)-!/'-h(1)) starting from the equation of mo- 
tion for tl^ni^)' 



dti 2m ^ ' 



dl'y(l-l')^|,(l>w(l')^H(l) = 0. (A15) 



Let us multiply eq. (|A15|) by (=F*/^)V'li(2) from the left, 
take its thermodynamic average, and compare the result 
with eq. (|A14[) . We then obtain the identity: 



dl'l7(l-l')(V|,(2)V|,(l')V'w(l')V'H(l)) 

±in J[Y.'^il,3)Gi2i3, 2) + Ei2(l, 3)G^(3, 2)] d3 . 

(A16) 



Setting 2 = 1 in eq. (|A16[) yields an expression of the inter- 
action energy in terms of the self-energy. An alternative 
expression is obtained by taking its complex conjugate. 
We thereby conclude that eq. (|A12p is indeed the inter- 
action energy of the system. 



2. Conservation laws 



Let us take complex conjugate of eq. l|A14p and use 
eqs. ([TIa)) . (|iga|) and dHaj). This yields 



OT2 



8 n'^xi'^ 



Gi2(l,2)-y"[G«(l,3)Si2(3,2) 
Gi2(l,3)S^(3,2)Jd3 = 0. (A17) 



2m 



First, we subtract eq. (|A17p from eq. (|A14[) . set 2 = 1, and 
use the identity (j A4bp . We thereby obtain the number 
conservation law as 



dn{l) 
dti 



Vi-j(l) = 0, 



(A18) 



where n{l) and are the particle and current densities 
defined by 



n(l) = ±i;iGi2(l,l), 



j{l)^±h^Yl_ZlGi,{l,2) 
2m 



(Al9a) 



(A19b) 



respectively. 

We next operate =F*ft(Vi — V2)/2m to eqs. (|A14[) and 
(|A17|) . subtract the latter from the former, and set 2 = 1. 
We thereby find that the time evolution of the current 
density obeys 

7^j(l) + -Vie^(l) + ^Vi[/(l) = 1q(1), (A20) 
oil m m m 



where Q(l) is given by eq. 
defined by 



and tensor 6^(1) is 



= T^(Vi.-V2.)(Vi,-V2,)Gi2(l,2)|2^^ . (A21) 



4m 



Integrating eq. (jA20p over the whole space of the relevant 
system and using eq. (jA7p . we obtain the total momen- 
tum conservation law as 



_d 

dt 



-/j(l)dVi = - /^Vi;7(l)dVi. (A22) 
iJ J m 



We finally operate ^ih-^ and ^ih-^ to eqs. (|A14p 
and (jAl7j) . respectively. Adding the resulting equations 
and setting 2 = 1, we obtain 

' V,.i;(l) + t/(l)^=g.(l), (A23) 



dti 



where (3e(l) is given by eq. (|A13|) . and £k(1) and ^^(l) 
are defined by 



fK(l) =±^Vi-V2Gi2(l,2) 

2m 



(A24) 



We are now ready to prove that the conservation laws 
are automatically satisfied in the ^-derivative approxi- 
mation. 



, (A25) 
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respectively. The quantity £k{^) denotes the kinetic en- 
ergy density. We next integrate eq. (|A23p over space and 
use eqs. ()A11| and (jAlSp . We thereby obtain the total- 
energy conservation law as 



_d_ 

dti 



£:K(i)dVi + (i?i„t(ti)) 



= -j3{l)-ViU{l)d\i. 



(A26) 

Thus, the conservation laws are automatically satisfied 
in the ^-derivative approximation as eqs. ()A18|) . (jA22p 
and ((A26l) . 



APPENDIX B: EXPRESSION OF V 

We here present expressions of the vertex function (|35p 
within the second-order skeleton expansion by using eqs. 
((26)) and ((27)) . This will help us to understand the struc- 
tures of r. 

First, Tij^jj/ is obtained from eq. 
It may be written compactly as 



with eq. 



r(^4.,(ll',22') = l7(l-2)%J,,,v(f3W['5(l,l')'5(2,2') 
±<5(1,2')<5(2,1')]. (Bl) 

Thus, only r^"^\;j^ = — r22^22 finite among T^^} 
(1'] 

Second, r--/^-^, is calculated from eq. ((77)) . It turns out 

that only F-^^-^ and F-^^^- are finite. For i^j, they are 
given by 

Fg^.(ll',22') = -^hV{l-2)V{l'-2')[G,,{l'^) 



xGj,(2',2)±Gj,(l',2)G,,(2',l)] 



(B2a) 



rS,(ii',22') 



^~ih F(l'-2)^(l-2')G,,(l',l)G,,(2',2) 

-f5(l',2)(5(l,2') jdij d3'T/(l'-3)F(l-3') 

X Gji(3, 3')Gij (3', 3) 
±5(l',2)F(l-2') / d3t^(l'-3)G,y(2',3)Gj,(3,l) 



±F(1'-2)5(1, 2') / d3 F(1-3)G,,(1', 3)G,,(3, 2) 



(B2b) 



respectively. Also, F^^ is given in terms of the above 
quantities as 



p(2) 
^ 11,11 


(11', 22') = 


-T%,{1%2'2) 


= -0{A 




h)T^^l,^{n\22') 




~h)e{t2 - 


QT%,{n' ,22') 




- t[)9{t2 - 


^2)r^i!2i(ll',22') 


-0{h 




h)T^^l,^{n' ,22') 



(B2c) 



In) 

These are the expressions of r^,-,^y for n = 1,2 in the 
space-time coordinates. 

We next write down F^^/ (11', 22') of the local ap- 
proximation which will be necessary later. In this case, 
Tiii J ji {11' , 22') acquires the structure of the uniform sys- 
tem. To be specific, we substitute eq. (H71) into eq. ((55)) 
and adopt the local approximation. We then find order 
by order that Fi^/ jj/ (11', 22') can be expanded as 

Fii' j j' (11', 22') 

d^pidei f d^padez f dVl^ »(pi+.r-i-ei+ti)/;i 



xe 



{2TTh)'^ J (27rft)4 J (27rft)4 

i[{pi--r[-ei-t\)-{p-2--r-2-e2-t2) + {p2+-r'^-e2+t'2)\/h 



xT,i' [piei, P2e2; quj, rt) , 



(B3) 



where Pj± =Pj±q/2 and ej±=ej±uj/2 in this expression, 
and rt denotes a local space-time point around 1, 1', 2 
and 2'. Let us write down T^^^ ^j,(pi,p2, q) explicitly for 

n = l,2 withp = pe and q = quj. First, T']^^\-^{pi,p2,q) is 
obtained from eq. ((Bip as 



(B4) 



rn4i(Pi>-P2, q) ^Vq± Vp^-p^ . 



Equations ((B2a)) and ((B2b)) for i ^ j arc transformed 
similarly into 



/d^q' 
j^^^Gji{pi-q')Gji{p2+q') 

(B5a) 



(27rft)4 

XVq'+q/2(Vq'_q/2±Fp2-pi+q'+q/2) 



r-Jj»(pi,P2,g) = 

x(K,'-q/2±V'pi- 

Also, eq. <\B2c\ yields 



-^^G,,{p^ + q')G,j{p2 + q') 



.(2) 



)(^q'+g/2±^pi-P2) 



rl2!i2(pi£i.P2£'2.g^) 

{ei+-e'^){e2+--e'.2) 

^?l2l{Pl£'l^P2£'2^q^) , ^21.2l{Pl£'l^P2£'2^q^) 



(B5b) 



/•de'i 


/•d4 " 


/ 2^ j 


27r _ 



(£i+-e;)(e2--£2) 

r2l!l2(Pie'l:P2£2,Q^) 



(£i_-e;)(e2+-£2) 
with £_,-!- =ej±?'0-|- in this expression. 



(£i_-£'i)(e2--e2) 

(B5c) 



APPENDIX C: CONSERVATION LAWS IN THE 
WIGNER REPRESENTATION 



Equations ((A18|) . ((A20() and ((A23p have a fundamental 
importance of describing the flows of particle, momentum 
and energy. We here transform these differential conser- 
vation laws into the Wigner representation within the 
first-order gradient expansion. 
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Substituting eqs. (|47)) and (|51ap into eq. (|A19p . we 
can write n(rt) and j{rt) alternatively in terms of A and 
(p. Those expressions are formally exact and satisfy eq. 
(|A18|) . The density n{rt) and the local velocity v{rt) = 
j{rt)/n{rt) are now given by 



n{rt) = h 



cL'^p de 
(27rft)4 



A{pe, rt)(t){pe, rt) , 



(Cla) 



virt) = ^ / ^Aipe, rmpe, rt) , (Clb) 
n[rt) J (zTThp m 



respectively. The particle conservation law (|A18|1 then 
reads 



dn 
9t 



V{nv) = 0. 



(C2) 



We next consider eq. (jA20[) for the momentum flow. 
Here it is desirable to express Q(l) on the right-hand side 
as a divergence. To carry this out within the first-order 
gradient expansion, we use eq. (jAlOP and its complex 
conjugate in eq. (jA8[) . We then obtain an alternative 
expression of Q{1) as 



Q(l) = 



with 1' = r'lti in this expression. We further write 
{i^\,{l)^|;\^{l')^n{l')i^n{l)) = P2{ri-r[,'^M) and 
expand the argument ni^l from ri up to the first 
order.'Sill We thereby obtain 



Q(i) = -Vin^(i), 



(C3) 



where tensor H is defined by 



n.^(i) 



1 /■ ,3-- dV{f) 
-— I d rr,- — 'xvn-iy 

2 df 



(V'|,(i+)V;,(i-)V^«(i-) 



x^w(l+)) , 



(C4) 



with l± = {ri±f/2,ti). Let us substitute eq. (jC3|) into 
eq. (|A20p . We then arrive at the differential momentum 
conservation law as 

^i{rt) + -Ve{rt) = -!^VU{rt) , (C5) 
Ot m m 



where tensor Qij is defined by 



e,,{rt)^efArt)+UYJrt) 



(C6) 



The quantity O^, given by eq. (|MT|) . may be written 
alternatively in terms of A and / as 



= mn{rt)v,{rt)vj{rt) + Uf^ (rt) , 



where 

n,^ (rt) ^ h J ^Mpe, rmpe, rt) , (C8) 

with p = p — mv . This {rt) denotes kinetic part of the 
momentum flux density tensor in the coordinate system 
moving with the local velocity v{rt). Using eqs. (jCip . 
((C2|) . ((C6|l and ([CT]) . we can transform eq. ((C5)) further 
into 



dv ^ 1 ^ VU 

— + v-Vv^ vn = 

ot mn m 

where tensor Ily is defined by 

Ii,,{rt)^Iif^{rt)+Ul{rt), 



(C9) 



(ClO) 



with and IlJ^ given by eqs. (|C8|) and (jC4[) . respec- 
tively. 

It remains to evaluate eq. \GA\ within the first-order 
gradient expansion of the ^-derivative approximation. 
This is essentially the calculation of the two-particle cor- 
relation function /C defined by eq. pO)) . Since li(Art) 



is operated by V in eq. (|C5|) . we can adopt the local 
approximation for this purpose. Thus, the procedure to 
obtain ICu'jj' (11', 22') is exactly the same as that for the 
uniform system. To be specific, let us substitute eqs. (|47|) 
and (jB3|) into eq. ([36|) . We then find that K. can also be 
expanded as 

^vi'.j]' (11', 22') 

d^Pldei f d^P2d£2 f ^^QAuJ r(p^_.r^-e^-t^)/h 

(27rft)4 J (27r;i)4 J (27rfi)4 

xlCu'.jj'{piei,p2e2;quj,rt) . (Cll) 
The quantity ICw jj'{pe,p'e';qLLj,rt) satisfies 



IC = ±(iTHiGGr) GG, 



(C12) 



where every quantity should be regarded now as a matrix 
in terms of pe-p'e' instead of ll'-22' in eq. ([M)) . with 
integration J d^pde/(2TTh)* over every internal variable 

pe implied. Indeed, 1 and GG are now defined by 

ii)u'.n'ipe.p's') = S,,6,,,,i2nh)^Sip-p')5{e-e') , 

(ClSa) 

{GG)u',jj'ipe,p'e';quj,rt) 

= G,j>{p.e.,rt)Gj,,{p+e+,rt){2nh)^S{p^p')S{e~e') , 

(C13b) 

respectively, with p± = p ± q/2 and e±=e ± u!/2 in this 
expression. Equation (jC12p enables us to calculate two- 
particle correlation functions for a given vertex function 
r, which in turn is specified completely for a given <I>. 



We now express 11]^- of eq. (|C4p in terms of the solution 
(C7) of eq. (|C12p . We first rewrite eq. (IC4| by using /C of eq. 
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([30)1 . We then substitute eqs. ©, (gT]) and (jClip into 
eq. (|C4p . remove by using e*'^ ''/^ = —ih-^e^'^'^/^^, 
perform partial integrations over q, and carry out the 
integration over f. We thereby arrive at an alternative 
expression of tensor 11^. Substituting it as well as eq. 



C8)) into eq. (jClOp . we obtain 



d^p'de' 
{2TTh)^ J (27rft)4 J {2Trhy 



x{Vq6ij+- 

q aq 

+ {2TTh)^S{q)S{iu)Gi2{pe,rt)Gi2{p'e',rt)] , (C14) 



where p — p^mv, and 6*12 and /C are given by eqs. (|51ap 
and (jC12p . respectively. We observe clearly that H is a 
symmetric tensor. 

We finally consider the differential energy conservation 
law. Equation (jA23p is not suitable for this purpose, 
however, because it is not written explicitly in terms of 
the local energy density. We hence start from the energy 
density defined by 

+ 1 J d'r[ y(ri-ri)(^t^(l)^|,(l')^«(l')^«(l)) , 

(C15) 

with t[ = ti in this expression. Let us differentiate eq. 
(ICISP with respect to time, eliminate time derivatives of 
the field operators by using eq. (jAlSp . and carry out the 
first-order gradient expansion for the interaction term. 
These standard proceduresSiil lead to the differential en- 
ergy conservation law: 



d£{rt) 
dt 



(C16) 



where denotes the energy flux density defined by 
3,{rt) = ±-^(V-V')V-V'Gi2(ri,r'i) 

+ ]- [ dfVif){^lir^t)jir+t)^Hir-t)) 



1 f^3--dV{r) 



■ Wn{r-t)j{r+t)i:n{r-t)) 
+ {^Ur+t)j{r_t)i;n{r+t))] , (C17) 



with 



j{rt) 



2mi 



(V-V')V4(r't)V'«(r<)L 



(C18) 



Since £ and in eq. (jC16p are operated by dt and V, 
respectively, eqs. (jClSp and (|C17p should be evaluated 
with the local approximation in the first-order gradient 



expansion. This can be carried out with the same proce- 
dures as those of deriving eq. (|C7p and the second term 
in eq. (|C14p from eqs. (|A2ip and (|C4|) . respectively. We 
finally obtain 



£ — —rnnv^ + £ , 



je — -mnv V + £v + Ilv + jq . 



(C19a) 



(C19b) 



Here tensor H is given by eq. (jC14p . and £ and jq are 
defined by 

£{rt) = h \ -JL^^A{ps,rt)^{ps,rt) 



2T:hY 2m 

+Y.i2{pe, rt)G^{pe, rt)] , (C20a) 



3Q(rt) = h I -^^pAipe, rt)(l){pe, rt) 



+ 



{2TTh)^ 2m?-' 
(iKf- f d^qdoj f d^pde f d^p'de' 



2T:hY J {2'KhY J (27rft) 



X - (2pV, + qSJ^^) [/Ci2.i2(pe, p'e'; qcj, rt) 
m\ q aq J 

+ {'2TTh)^d{q)d{u)Gi2{pe,rt)Gi2{p'£',rt)], (C20b) 

respectively, with p = p — mv. Use has been made of eq. 
(IA16P to derive the interaction term in eq. (|C20ap . The 
quantity £ denotes the energy density in the reference 
frame moving with the local velocity v, whereas jq is 
the heat-flux density. 

Let us substitute eq. (|C19P into eq. (|C16p and trans- 
form it with eqs. (|C2p and (|C9p . We thereby arrive at 
an alternative expression of the differential energy con- 
servation law as 

^ + V-{£v+jQ) + J2u,,^^=0. (C21) 



Equations ^ and (IC2T]) with eqs. dUT]), (ICMl) 

and (jC20p completely describe the flows of particle, mo- 
mentum and energy, respectively. 



APPENDIX D: CALCULATION OF $ IN 
EQUILIBRIUM 

Using the zero-temperature time-ordered Goldstone 
techniquCf^ Carneiro and Pethick^ performed a third- 
order calculation of the functional $ for a uniform Fermi 
system. They thereby found a singular or on-energy- 
shell contribution to $; see eq. (29) and §1V of their 
paper. It is this on-energy-shell contribution that brings 
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FIG. 5: Third-order diagram for 



the difference between eq. ((69|) and the Carneiro-Pethick 
expression. As they mentioned exphcitly in §IV.A, how- 
ever, there may be ambiguity in the Goldstone technique 
on how to regularize the energy denominators. We here 
calculate the same contribution to $ with the finite- 
temperature Matsubara formalism, choosing the chem- 
ical potential as an independent variable instead of the 
total particle number. The Matsubara formalism has def- 
inite advantages over the Goldstone technique that (i) it 
can describe finite temperatures and (ii) no additional 
regularization procedure is required. 

Let us express the Matsubara Green's function in the 
Lehmann representation as 



G{ri,r2,Zn 



de A{ri,r2,s) 
27r Zr, — s 



(Dl) 



where z„ = inirk^Ti for bosons and Zn = {2n+l)T:kBTi 
for fermions. We also introduce the bare vertex: 

(rir^lylnra) = y(n -^2) [5(r'i -ri)5(r^ -ra) 

±S{r[-r2)S{r'2-r,)] , (D2) 

which is expressed by a square in the Feynman dia- 
gram following the convention of Abrikosov, Gor'kov and 
Dzyaloshinski.'^^ With eq. (jD2p . all the third-order con- 
tributions to $ can be expressed by a single diagram of 
Fig. m Using eqs. (|D1|) and (ID2|) . we can write down the 
corresponding analytic expression compactly as 



Here X^^^ is defined by 

, £2, ^3, ^4, £5,^6) 

1 '^ni+n2,n3+n4'^ni+n2,n5+n6 
fl4 



(D3) 



„f::^8 (21 -£i) (2:2 -£2) (23 -£3) (2:4 -£4) 



(z5-£5)(26-£6) ' 



(D4) 



with /3 = \/k^T and z^ = z^-. The other factor K!^'^^ 
is analytic in £, and irrelevant for the present purpose. 



Indeed, it is given explicitly as 

, £2, £3, £4, £5, £6 J 

X (r;r^|y|rir2)A(ri, r^, £5)A(r5, r^, £3)A(r3, r^, £1) 
X [=F4>l(r2, r4, --£2)A(r4, rg, -£4)A(r6, r'^, -£5) 

+A(r2, r^, £6)A(r6, r^, £4)A(r4, r^, £2)] , (D5) 

where the first (second) term in the square bracket cor- 
responds to the ring and particle-hole (particle-particle) 
contribution. 

The whole issue here is whether the above Z^^-* contains 
the on-energy-shell contribution. However, the expres- 
sion (jD4p already tells us the absence of the on-energy- 
shell contribution for the Fermi system considered by 
Carneiro and Pethicki^l because Zj is pure imaginary and 
£j is real, the fraction is regular at any finite tempera- 
ture, even for £i-|-£2 = £3 + £4 = £5 + £6. We shall confirm 
this fact further below. 

We now carry out the summations over rij one by one 
with the standard technique of transforming them into 
contour integrals with^ 



1 



e/3^ =F 1 ■ 

First, those over na and 71.4 in eq. (|D4[) yield 

0304-(l±03)(l±<^4) 



(D6) 



-y 

^«tr. (^3-^3)(^4-£4) 
= ^^(^^(21,22, £3, £4) , 



Zi-|-Z2-£3-£4 



(D7) 



with 



In obtaining the result, we have used 



0(zi + Z2 — £4) = ^(-£4) = t(1±04) and l±(/)3±</)4 = 
(l±(/)3)(l±04)-(/)3(/)4. Note that J^'^'^ is regular for the 
imaginary arguments z\ and Z2. The summations over 
715 and ng can be performed similarly. 

Before proceeding directly to the summation over n2 in 
eq. (jD4p . it is useful to consider a couple of summations 
connected with eq. (|D7p . The first one 5*^^^ is defined 
and transformed as follows: 

yZXi Z2, £3, £4) 

5^^^ (01, £2, £3, £4) = -^2^ : — : 



02 -£2 



± 



(/)2(l±(^3)(l±04)T(l±'^2)(/'304 



^l+£2 — £3 — £4 



(D8) 



where we have used l±03±(/)4 = 0304[e'^^^^"'"^*'^ — 1] and 
[e'^(=3+'=-*)-l](/)(£3-H£4-zi) = ±l. The second one X^^) jg 
given by 



£2, £3, £4) 



-y 



5(2) I 



^l, £2, £3, £4) 



/3 



Zl-£l 



(l±(/.l)(l±02) 



^3<?^4 



62(1±<^3)(1±<^4) 



£l+£2-£3-£4 



(D9a) 
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where use has been made of the identity [4>2 ( l±i/>3 ) ( Ui/'4 )t 
(l±(/)2)^3'A4]</'(e3+£4-e2) = (l±02)^3'^4- Equation (|D9a[) 
is what we encounter in the second-order calculation for 
Note that X^^^ is analytic even when ei+£2 = £3+£4- 
Thus, we may write I^^-* alternatively as 

X(^)(£i,e2,e3,e4) 

^p (l±0l)(l±02)03^4 p (^l02(l±03)(l±04) pg^^ 



ei+£2 — £3^£4 



£1+62 — £3 — £4 



with P denoting the principal value. 

Now, after those over 113, n^i, and tiq given by eq. 
(|D7p . the summation over 71,2 in eq. (|D4[) is performed as 
follows: 



3) 



, £2, £3: £4, £51 Sej 
-(3^ Z2~e2 

= {5(2)(zi,e2,£3,e4)[</'5'/'6-(l±</'5)(l±06)] 
^5(2)(2l,£2,e5,e6)[03'^4-(l±03)(l±04)]} 

1 



£3 +£4 — £5 



(DIO) 



This iS'^'' is clearly analytic at £3 +£4 = £5 +£6. Finally, 
eq. (jD4p is transformed as 



, £2, £3, £4, £5: £6j 
_ 1_ 5('^)(zi,£2,£3,£4,£5,£6) 

~ zi-ei 

Til 

= {X(2)(£^,e2,e3,£4)[^50g„(l±05)(l±0g)] 
^X(2)(£i,£2,£5,£6)[</>304-(l±03)(l±04)]} 



£3 + £4 — £5 — £6 



(Dll) 



Remember that T^^^ is analytic as eq. (|D9aP . In addition, 
has no singularity at £3 + £4 — £5 +£6. We hence 
conclude that T^^^ is analytic, i.e., there is no on-energy- 
shell term in Z*^^^. This is clearly a general feature of $, 
as may be realized most easily from its expression with 
respect to the Matsubara frequency such as eq. (|D4p 



From eqs. (jPSP and (|D9bp . we obtain 



(5l(2)(£i,£2,£3,£4) (2)/ X 
— = TRc5^''^(£i+,£2,£3,£4) : 



(D12) 



with £1+ =£i-|-zO+. The same relation holds between l'^^ 
and S^'^\ Using them, we obtain eqs. (3)-(5) of ref. H^- 
It should be noted that the same expression as eq. (5) of 
ref. had been presented by Fulde and Wagner'*- for a 
Bose system without any detailed derivations. 
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o 
o 



In a previous paper I have proposed a principle of max- 
imum entropy for nonequilibrium steady statesii How- 
ever, the argument based on the H theorem in Intro- 
duction is not appropriate, because there is net outflow 
of entropy in nonequilibrium steady states. Thus, the 
principle is to be regarded as a pure conjecture. 

It should also be noted that the principle is identical in 
the dilute classical limit to the information theory devel- 
oped by Jou and co-workers.-^ Kim and Hayakawai^ tested 
the information theory for dilute classical hard-core and 
Maxwell molecules with heat conduction by comparing 
its predictions with those from the steady-state Boltz- 
mann equation. They found that the two approaches 
can yield identical results only up to the first-order devi- 
ations from the local equilibrium generally, and discrep- 
ancies emerge at the second order for all the physical 
quantities except entropy. This implies that the above 
principle cannot be valid beyond the first order. It is pos- 
sible, however, that it generally holds within this lowest 



order, as exemplified by Kim and Hayakawa'^ and also 
by ref. 4 on Rayleigh-Benard convection. Another sup- 
port to this conjecture is that the nonequilibrium statis- 
tical operator of Zubarev corresponds to an extremum 
of the information entropyi^ Note finally that, unlike the 
linear response theory which is an expansion from the 
global equilibrium, the first-order expansion from the lo- 
cal equilibrium can describe a wide range of nonequi- 
librium nonlinear phenomena, as may be realized by the 
fact that the Navier-Stokes equations belong to this order 
of approximation. We certainly need further investiga- 
tions to clarify the applicable range of the conjecture. 



Acknowledgments 



I am grateful to H. Hayakawa for critical discussions. 



Soc. Jpn. 75 (2006) 114005; 



1 T. Kita: J. Pliys. 
cond-mat/0611270 . 

^ D. Jou, J. Casas- Vazquez and G. Lebon; Extended Irre- 
versible Thermodynamics (Springer- Verlag, Berlin, 1993). 

^ Kim H.-D. and H. Hayayawa: J. Phys. Soc. Jpn. 72 (2003) 
2473. 

T. Kita; J. Phys. Soc. Jpn. 75 (2006) 124005; 



cond-mat/0611270 

D. N. Zubarev: Nonequilibrium Statistical Thermodynamics 
(Consultants Bureau, New York, 1974). 
S. Chapman and T. G. Cowhng: The Mathematical Theory 
of Non-uniform Gases (Cambridge University Press, Cam- 
bridge, 1990). 



